Colloidal electroconvection in a thin horizontal cell: I. microscopic cooperative 

patterns at low voltage 
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Applying an electric field to an aqueous colloidal dispersion establishes a complex interplay of 
forces among the highly mobile simple ions, the more highly charged but less mobile colloidal spheres, 
and the surrounding water. This interplay can induce a wide variety of visually striking dynamical 
instabilities, even when the applied field is constant. This Article reports on the highly organized 
patterns that emerge when electrohydrodynamic forces compete with gravity in thin layers of charge- 
stabilized colloidal spheres subjected to low voltages between parallel plate electrodes. Depending 
on the conditions, these spheres can form into levitating clusters with morphologies ranging from 
tumbling clouds, to toroidal vortex rings, to writhing labyrinths. 



INTRODUCTION 

Electric fields exert forces on charge-stabilized colloidal 
particles both directly through their coupling to the par- 
ticles' charges and also indirectly through their influence 
on the surrounding electrolyte. Neighboring particles 
also interact electrostatically and hydrodynamically with 
each other as they move. The resulting cooperative mo- 
tions in many-body suspensions can be quite complex, 
and our understanding of such electrokinetic phenomena 
remains incomplete despite more than a century of study 


The tendency of colloidal particles to form oriented 
chains in oscillatory (AC) electric fields was noted nearly 
a century ago and was correctly ascribed to field- 
induced dipolar interactions. More recently, AC fields 
also have been found to organize colloidal spheres into cir- 
culating chevron bands oblique to the fieldjp,0,0,0,0, 
glassy bands perpendicular to the field |8|, and even 
highly ordered colloidal crystals H,EiE3- Still other pat- 
terns form when electric fields are used to drive charged 
colloid against electrodes' surfaces, with two-dimensional 
fractal aggreg ates appeared in monodisperse suspensions 
[miT^ITa. 14| and planar superlattices forming in binary 
colloids [15 1 . More recently, a colloidal model system has 
been introduced whose phase diagram can be tuned with 
salt concentration and AC electric field (lif. 

Even constant (DC) electric fields give rise to complex 
dynamics. For example, charged spheres driven onto elec- 
trodes' surfaces by DC fields self-organize into epitaxial 
colloidal crystals 13 E 113 DC fields also have 
been observed to induce complex labyrinthine patterns 
in nonaqueous nanocolloids |21|. 

Most of these cooperative electrokinetic effects have 
been explained on the basis of electrohydrodynamic cou- 



pling among the spheres mediated by field-induced fluxes 
of ions. Others, such as the electrokinetic compression of 
colloidal crystals |9( have so far resisted explanation. 

In this paper, we report the remarkably diverse set of 
patterns that charged colloidal particles can form as sedi- 
menting under gravity in a vertical DC electric field. Pre- 
liminary results have been published in Ref. [2^ . Beau- 
tiful, highly organized dynamic patterns with spatial pe- 
riods between 20 and 200 micrometers form at biases 
near the threshold for electrolysis. These are supplanted 
by macroscopic (spatial period ~ 0.5 — 2 mm) patterns 
at higher biases. We describe how such patterns form 
through the interplay of electrohydrodynamic coupling 
driven by reaction-diffusion of ions and the uniform body 
force provided by gravity. We further show that the mi- 
croscopic patterns result from many-particle interactions 
rather than an underlying convective instability of the 
electrolyte. 



EXPERIMENTAL SYSTEM 

The experimental apparatus is shown schematically in 
Fig. n] Our samples consist of aqueous dispersions of 
silica spheres 3.0 fj,m in diameter (Bangs Laboratories, 
Lot Number 4181) confined to slit pore between a glass 
microscope slide and a coverslip. The inner glass sur- 
faces were coated with 10 nm thick gold electrodes on 
10 nm thick titanium or chromium wetting layers before 
assembly. While still optically thin, the electrodes have 
a resistivity of less than 50 ohms per square and allow us 
to apply uniform vertical electric fields to the confined 
suspension. Conductive epoxy is used to make contact 
to the electrodes once the cell is assembled. In most 
cases, glass spacers were used to set the electrode sepa- 
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FIG. 1: Parallel-plate electrochemical cell with semitranspar- 
ent electrodes on the inner glass surfaces. The fluxes of ions 
are shown schematically together with the electrode reactions 
responsible for their evolution. A single colloidal vortex ring 
is shown schematically, levitated to height h in the inter- 
electrode gap H . The arrows indicate the vortex's observed 
circulation direction. 



ration, H , in the range 100 fim < H < 400 /jm across the 
4 x 1.5 cm 2 observation area. Glass access tubes bonded 
to holes drilled through the slides were left open to air so 
that the suspension equilibrates to a pH of roughly 5.5 
through the dissolution of carbon dioxide. Under these 
conditions, the silica spheres acquire a surface charge 
density of roug hly -0.4 mC/m 2 [2J. As the density of 
colloidal silica is roughly 2 g/cm 3 , the 3.0 /mi diameter 
spheres sediment rapidly onto the the lower electrode. 
Consequently, we use the areal coverage <f> to report the 
concentration of particles, with (j> = 100% corresponding 
to a close-packed monolayer and 4> = 200% correspond- 
ing to a bilayer. Colloidal patterns are imaged with CCD 
(charge-coupled device) camera through an inverted op- 
tical microscope. The images and videos reported in this 
paper were obtained looking upward through the lower 
electrode, in the direction opposite to gravity. We define 
the applied DC voltage to be positive if the upper elec- 
trode is positive so that the electric field is directed down- 
ward, but the electrostatic forces on negatively charged 
particles point upward. 

The electrochemical processes in a parallel-plate wa- 
ter electrolysis cell are described in most standard refer- 
ences, such as Refs. [24| and [2^, and are summarized in 
Fig. n Decomposing a water molecule requires an input 
of 1.23 eV, independent of the salt concentration or the 
pH value of the solution. In practice, however, the redox 
potential is at least 0.6 V higher due to a voltage drop 
across the electrode-electrolyte interface known as the 
overvoltage (overpotential). Once the threshold for hy- 
drolysis is reached, the electrodes rapidly become coated 
with insulating layers of H2 and O2 . These barriers soon 
stop the electrode reactions unless the evolved gases form 
bubbles that escape the electrodes. A substantial nucle- 
ation barrier impedes the evolution of H2 and O2 bubbles 
large enough to expand and escape the electrodes. Con- 
sequently, a bias as large as 2.5 V to 3.0 V typically is re- 
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FIG. 2: Current density j(V) as a function of applied volt- 
age in samples of deionized water in three cells of thickness 
H = 90, 200 and 370 (jm. A single cycle of increasing and de- 
creasing bias is shown for each sample. Values were recorded 
after 130 sec relaxation at each voltage. The intersection of 
the extrapolated linear trends yields the threshold for hydrol- 
ysis. 



quired to maintain hydrolysis in quasi-steady-state. The 
range of threshold biases is due to inevitable variations 
in the surface condition of the electrodes. 

Our sample cells' vacuum-deposited thin film gold elec- 
trodes yield threshold biases between 2.6 and 3.2 V. Vis- 
ible bubbles begin to appear between a minute and an 
hour after the bias is applied, depending on the volt- 
age. The onset of bubble evolution sets our observation 
time, and decreases with increasing voltage or salt con- 
centration. Although the redox potential for gold in wa- 
ter is only 1.5 V, the overpotential prevents the anode 
from etching, so that a negligible concentration of gold 
ions, Au 3+ , is released into solution. Indeed, no visi- 
ble or electrochemical change in the electrodes is seen 
even after several hours operation at a steady bias. Even 
in the absence of added salt, however, the electrolyte 
in our samples has total ionic concentrations of roughly 
n HCO~ ~ n H 3 0+ ~ 3 A«M and jioh- ~ 3 nM. A 0.1 /im 
thick Debye-Hiickel screening layer forms at each elec- 
trode at these concentrations, beyond which the field 
would not penetrate at equilibrium. During hydroly- 
sis, however, the electrodes act like sources and sinks 
of H3O" 1 " and OH - ions whose currents establish both 
an electric field and a pH gradient through the cell, with 
the lower half volume becoming acidic and the upper half 
basic under positive bias. 

To quantitatively characterize hydrolysis in our cell, we 
measured the current density j(V,t) under the simplest 
conditions: the parallel- plate cell filled with deionized 
water without colloid. We performed sweeping voltam- 
metry by increasing the applied bias in discrete steps 



3 




V(V) 

FIG. 3: Current density as a function of applied voltage at 
different salt concentrations. The samples are aqueous solu- 
tions of NaCl without colloid in a cell H = 200 )im thick. The 
current was measured after 130 sec relaxation at each step in 
an increasing sequence of applied biases. 



of 0.2 V from V to 3.6 V and then back down to V, 
waiting between 2 to 5 minutes at each voltage level. Fig- 
ure [21 shows the voltage dependence of the steady-state 
current in cells of different inter-electrode separations, H. 
Each curve consists of two linear parts, one at low bias 
in which little current flows, and another ohmic regime 
at higher biases. Their intersection yields the threshold 
for water decomposition |25|. Below the decomposition 
voltage, the small but nonvanishing current is due to the 
continuous dissolution of H2 and O2 molecules. This is 
known as the residual current and also increases with the 
applied voltage. 

The higher conductivity of the high-bias regime is fed 
by the additional ions from the decomposition of water. 
Adding NaCl increases the conductivity but does not af- 
fect the threshold, as shown in Fig. [3] For this reason, 
NaCl is referred to as indifferent electrolyte. 

From the slopes of the data in Fig. [21 we estimate 
the equilibrium conductivity of our electrolyte to be 
a ss 4 x 10~ 4 n _1 m . Adding colloidal spheres to 
the deionized water increases the ionic concentration and 
thus the conductivity. Adding roughly <j> = 100% of silica 
or polystyrene sulfate (PS) colloid increases the conduc- 
tivity to 7 x 10~ 4 and 20 x 10~ 4 _1 m _1 , respectively. 
From these values, we estimated the pH for these systems 
to be 5.0, 4.7 and 4.2, respectively, with the negatively 
charged colloidal spheres acting as weak acids. 



COLLOIDAL ELECTROKINETIC PATTERN 
FORMATION 

Interfacial crystals 

The field-driven fluxes of ions and charged particles are 
coupled to flows in the surrounding electrolyte and can 
give rise to electrohydrodynamic instabilities. The most 
commonly observed manifestation is the formation of col- 
loidal crystals on the electrodes 0]. Charged colloidal 
spheres are driven toward the oppositely charged elec- 
trode, and might be expected to diffuse independently 
on the surface. Instead, the particles rapidly coalesce 
into close-packed epitaxial crystals [l^. Th e strong and 
long-ranged interparticle attraction responsible for this 
crystallization results from an attractive hydrodynamic 
coupling induced by the spheres' distortion of the local 
fluxes of simple ions. 

Such electrohydrodynamically bound crystals also can 
be induced by an AC field 26]. At high enough fre- 
quencies, however, the dielectric spheres' induced dipole 
moments mediate an in-plane repulsion strong enough 
to and the interfacial crystals dissociate destabilize the 
interfacial crystals j2(|. 

Field-induced surface crystals have been observed in 
silica, PS, gold and even neutral PS colloidal spheres with 
a wide range of diameters, and on both indium tin ox- 
ide (ITO) and gold electrodes. The threshold voltage for 
crystallization can be very low if the spheres' diffusion is 
slow and if their inherent electrostatic and dipolar repul- 
sions are relatively weak. 

In our system, we have observed that the threshold 
voltage for field-induced crystallization depends on the 
particles' mobility and charge as shown in Fig. 0] but 
does not depend on the thickness of the cell and can be 
lower than the 1.23 V redox potential of water. Crystals 
were formed on the lower electrode with a negative DC 
bias. Larger spheres crystallize at lower voltages both be- 
cause they are less diffusive and also because the larger 
distortions they cause in the surrounding ionic fluxes give 
rise to stronger attractions. Nominally neutral hydroxyl- 
terminated polystyrene spheres reproducibly crystallize 
at the lowest voltages of the particles we studied, presum- 
ably because of their substantially weaker electrostatic 
repulsions. 

Colloidal vortex rings 

In this section, we examine some more complex phe- 
nomena that arise under a more restricted set of condi- 
tions, and which provide new insights into the interplay 
of electrohydrodynamic and other forces in macroionic 
systems. 

Applying a positive bias in our cell tends to drive neg- 
atively charged silica spheres upward into the bulk elec- 
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FIG. 4: Threshold voltage for interfacial crystallization for sil- 
ica, polystyrene sulfate (PS), and neutral polystyrene spheres 
over a range of particle diameters. 

trolyte. The field is uniform in our parallel plate geome- 
try, so that the upward electrophoretic force on spheres 
of constant charge therefore should be independent of 
height h in the cell. Under such conditions, spheres dense 
enough to sediment under gravity should remain near the 
lower electrode at low biases, and should rise directly to 
the upper electrode at higher biases where electrohydro- 
dynamics forces exceed gravity. 

This effect should be accentuated in charge-regulating 
silica spheres because the magnitude of their surface 
charge increases with pH, and therefore should increase 
with h. Indeed, individual isolated silica spheres ex- 
hibit hysteresis in their promotion to the upper electrode 
with increasing bias, and relaxation to the lower elec- 
trode with decreasing bias, consistent with these con- 
siderations. More concentrated dispersions, on the other 
hand, exhibit substantially more interesting collective be- 
havior. 

At biases between 2.2 and 4 V, which only slightly 
exceed the threshold for hydrolysis, the background elec- 
trolyte is stable against convective instabilities. Never- 
theless, initially quiescent collections of spheres develop 
several classes of striking many-body convective patterns 
that highlight a previously unsuspected role for cooper- 
ativity in charge-stabilized dispersions subjected to elec- 
tric fields. Unlike the electrohydrodynamically stabi- 
lized surface crystals described in previous reports, these 
quasi-steady-state patterns form in the bulk of the elec- 
trolyte and reflect the interplay of electrohydrodynamic 
forces and gravity. 

Figure shows the typical sequence of dynamic pat- 
terns that forms when an upward bias is applied to dis- 
persions of 3.0 /im diameter silica spheres in water. Pos- 
itively biasing the upper electrode by as much as 0.8 V 







| t_i l__mti__ im I - ■ ■ ■ — * 


t 


ftfA 



FIG. 5: Typical patterns formed in dispersions of 3.0 /jm 
diameter silica spheres at H = 200 ^m. (a) V = 0: Equilib- 
rium monolayer, <f> ~ 100%. (b) V = 2.6 V: Diffuse tumbling 
clouds, (c) V = 3.0 V: Circulating flower-like cluster levi- 
tated to h = 40 lira, (d) V = 3.8 V: Compact circulating 
cluster, also at h — 40 /im. 



causes no discernible out-of-plane motion. Abruptly ap- 
plying a larger bias in the range 0.8 V < V < 2.5 V 
causes a transient in which spheres jump off the lower 
surface and then settle back to the bottom. This tran- 
sient reflects the establishment of static ionic gradients 
within the water that screen out the field. 

The spheres' collective behavior passes through dis- 
tinct regimes as they are driven out of equilibrium. Just 
above threshold, the monolayer of spheres breaks into 
diffuse billowing clouds occupying the lower half of the 
sample cell, as shown in Fig. OJb) . Pushing the system 
further from equilibrium might be expected to yield in- 
creasingly chaotic behavior; but quite the opposite oc- 
curs. Increasing the bias beyond 2.6 V coalesces the itin- 
erant clouds into extraordinary flower-like clusters such 
as the example in Fig. [Sic), all floating with their mid- 
planes at h — 40 /im above the lower wall. Each cluster 
forms around a rapidly circulating toroid whose spheres 
travel downward along the inner surface and return up- 
ward along the outside, completing one cycle in a few sec- 
onds. The sense of this circulation is indicated schemat- 
ically in Fig. n] Because these clusters' structure and 
motion are reminiscent of hydrodynamic vortex rings, we 
will refer to them as colloidal vortex rings. 

Colloidal vortex rings formed at higher biases tend to 
consist of fewer clusters, each containing a larger num- 
ber of spheres at higher density. The dense yet vigor- 



ously circulating cluster in Fig. Etd) was observed at 
V = 3.8 V. Increasing the bias once clusters have al- 
ready formed causes each to shrink and its circulation 
to accelerate. Such over-driven clusters merge with their 
neighbors until they achieve the size appropriate to the 
final bias. The clusters' height, h, does not change ap- 
preciably with voltage. 

By around V — 4 V, some clusters' density increases 
enough that they become jammed and stop circulating. 
These compactified clusters still dissociate immediately 
once the bias is removed. 

At higher voltages, the electrolyte itself becomes un- 
stable against electroconvection. The resulting patterns 
span the electrochemical cell, with the colloid acting prin- 
cipally as passive tracers. We will discuss this distinct 
range of conditions elsewhere. 

Most often, a densely packed colloidal vortex ring is 
surrounded by a diffuse circulating corona that extends 
outward for tens of micrometers. Fig. [fj] shows alternate 
forms consisting of compact circulating toroidal clusters 
without coronas. Such bare colloidal vortex rings form 
most often at slightly lower voltages, 2.6 V versus 2.8 V, 
and can coexist with fully dressed colloidal vortex rings. 
These clusters also are levitated to h » 40 ^m in a system 
of thickness H — 200 fim. 




FIG. 6: Bare colloidal vortex ring with and without cen- 
tral core. 3.0 /im diameter silica spheres in H = 200 fim 
cell, (a) Compact toroidal cluster enclosing a domain of two- 
dimensional colloid crystal at V = 2.8 V. The toroid com- 
pletes one cycle of circulation in about 1.5 sec. (b) Compact 
toroidal cluster with empty core. 

Although colloidal vortex rin gs s omewhat resemble 
conventional laminar vortex rings |27j , they are driven by 
quite different mechanisms [2^ and also have other dis- 
tinctive features. Perhaps their most remarkable charac- 
teristic is the presence of close-packed colloidal crystals 
at their centers. The example in Fig. [BTal is particu- 
larly noteworthy because much of its core consists of a 
static crystalline monolayer. These close-packed spheres 
are not flocculated, and disperse immediately once the 
driving field is turned off. Consequently, the crystals' 
formation and stability suggests the existence of a stag- 
nation plane along each cluster's midplane, quite unlike 
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FIG. 7: Central crystal formed by compression in a disper- 
sion of 3.0 fim diameter silica spheres in a layer of thickness 
H = 200 fim. The sequence of photographs show a cluster 
shrinking as the bias increases from 2.4 V to 4.0 V. An AC 
component (0.7 V square wave at 10 Hz) is applied to further 
stabilize the crystal. The circulation period of the toroid is 
about 1.5 sec. (a): sec. (b): 0.5 sec. (c): 1.0 sec. (d): 
1.5 sec. (e): 2.0 sec. (f): 3.0 sec. 

the streaming flow within a purely hydrodynamic vortex 
ring. 

The stagnation plane must end at the edge of the cir- 
culating toroid, and indeed particles from the core some- 
times are swept into the circulating flow. While the crys- 
talline domains tend to be stable for the duration of an 
experimental run (up to ten minutes at V = 3 V), some 
are incorporated into the surrounding toroid after a few 
minutes, leaving hollow circulating rings such as that in 
Fig. Eft). 

In most cases, the crystalline monolayer core evolves 
from a layered fluid-like state as shown in Fig. As the 
ring's diameter decreases with increasing bias, the con- 
centration of spheres in the core increases to the freezing 
point, and a static crystal results. 

Individual colloidal vortex rings sometimes develop 
breathing-mode instabilities with periods of a few sec- 
onds, with one cycle of a typical structural oscillation 
appearing in Fig. [SJ The usual circulation continues 
throughout the oscillation. Remarkably, the cluster's 
crystalline core is destroyed and reformed with each cy- 
cle. Neighboring clusters' oscillations do not become 
phase locked, nor do they necessarily pulsate at the same 
frequency. Indeed, oscillating clusters can coexist with 
steadily circulating or compactified clusters. 

Within the range of biases for which flower-like col- 
loidal vortex rings reproducibly form (about 2.6 to 3.4 V), 
we also observe a variety of variant vortex ring structures 
that either coexist with or even replace the morphologies 
represented in Figs. [S] and Some of these are shown in 
Fig. El 

Figure EJa) shows colloidal vortex rings whose coronas 
have sharply defined boundaries, rather than the diffuse 
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FIG. 8: This sequence of images, separated by 0.5 sec, shows 
one period of a vortex ring's breathing mode oscillation. 
3.0 fim silica spheres in H — 200 /im cell at 2.8 V constant 
bias. The cluster's core remains crystalline throughout the 
cycle. 




diffusely bounded clusters that undergo breathing mode 
instabilities the corona in the diffuse morphology can un- 
dulate like a jellyfish. Figure |jjc) shows an extreme ex- 
ample of this circumferential instability. 

In all of these cases, the vigorously circulating vortex 
ring appears to be virtually decoupled from the central 
crystalline core. Figure Efd) shows a relatively rare case 
in which some of the particles in the central core are 
swept into the vortex ring's circulation. 

In roughly half of our experiments, we observed vor- 
tex ring s coexisting with two-dimensional surface crystals 
U3 Il8l | that had nucleated on the upper electrode. The 
spheres in these crystals are driven to the upper elec- 
trode by the transient ionic current burst when the bias 
is initially applied. Particles that reach the upper elec- 
trode during these transients tend to stay there and sub- 
sequently form crystalline clusters thr oug h mechanisms 
that have been previously reported 0,0]- Indeed, such 
surface crystals begin to form at biases below the onset 
for vortex ring formation. Above 2.6 V, clusters hover- 
ing in the bulk coexist with crystals trapped on the upper 
surface. Decreasing the bias to 1.4 V, which is the thresh- 
old for forming crystals on the upper electrode, causes all 
particles in crystals to sink to the bottom. 



Worm-like colloidal vortices 

Not all colloidal vortices adopt a ring-like configura- 
tion, as the examples in Fig. 1101 show. Lower concen- 
tration suspensions tend to favor elongated vortex rings, 
such as those in Fig. llOf a). that can extend into linear 
worm-like structures such as those in Fig. HUf b). These 
clusters continue to circulate, but as linear rolls rather 
than vortex rings. As a result, worm-like colloidal vor- 
tices interact more strongly that colloidal vortex rings, 
and can organize themselves into large-scale labyrinthine 
patterns, as in Fig. HUf cl. As for the other patterns we 
have described, worm-like colloidal vortices are levitated 
to roughly h — 50 in an H = 200 /im cell. 



Sedimented interfacial patterns 



FIG. 9: Variations on the theme of colloidal vortex rings. 
3.0 /im diameter silica spheres in H = 200 /im cells, (a) 
3.0 V: Colloidal vortex rings with sharply defined coronas. 

(b) 2.6 V: Sharply bounded vortex ring without central core. 

(c) 2.6 V: Central crystal stabilized by a 0.1 V square wave at 
1.0 Hz develops scalloped boundary, (d) 3.0 V: Rare variant 
in which the central core also is convective, cj> ~ 150%. 



coronas evident in Fig. [5Jc) . These clusters also some- 
times form without a crystalline core as can be seen in 
Fig. infb), and the two types can coexist. Unlike the 



Although the patterns described in the previous Sec- 
tion form reproducibly, they occasionally are replaced 
by qualitatively distinct labyrinthine and island-like pat- 
terns such as those shown in Fig. ^] Unlike the vig- 
orously circulating free-floating colloidal vortices, these 
slow-moving clusters form near the bottom electrode un- 
der the same conditions. These patterns also are distinct 
from the previously reported interfacial crystals, which 
form on the upper electrode under positive bias. 

Once sedimented interfacial clusters nucleate, the 
choice between forming islands, as in Fig. Illf a - ) and 




FIG. 10: Worm-like colloidal vortices, (a) 3.4 V: elliptical 
colloidal vortex rings, both with and without crystalline cores, 
coexisting with worm-like colloidal vortices, (b) 3.4 V: worm- 
like vortices extending from vortex rings, (c) 2.6 V: labyrinth 
of worm-like vortices. 



Illf bL or labyrinths, such as Fig. Illf c) and Illf cO. de- 
pends principally on the concentration of particles. Two- 
dimensional islands form at surface coverages of 10% < 
4> <■ 60%, and two-dimensional labyrinths at 60% < <\> < 
80%. Three-dimensional labyrinths form at higher con- 
centrations. Islands and labyrinths can coexist at inter- 
mediate sphere concentrations. 

The spheres in islands and labyrinths usually flow flu- 
idly as the patterns form, in a manner that resembles 
the circulation in bulk colloidal vortices. Unlike bulk 
colloidal vortices, whose circulation retains its sense over 
time, labyrinthine domains sometimes flip circulation di- 
rection as they coarsen. Unlike Rayleigh-Benard convec- 
tion or bulk electroconvection, moreover, neighboring do- 
mains' circulation directions generally are not correlated. 
Larger labyrinthine domains, however, develop double- 
roll structures, with particles rising along the edges and 
sinking along the center line. 

After a few minutes of circulation, islands and 





FIG. 11: Sedimented interfacial patterns, (a) Islands of 
3.0 /im diameter silica spheres in H = 90 /im cell at 2.6 V. 
(b) The same sample, after increasing the bias to 2.8 V. (c) 
Labyrinths of 3.0 /im diameter silica spheres in H — 200 fim 
cell at 2.6 V. (d) The same sample, after increasing the bias 
to 3 V. 



labyrinths often become jammed. Increasing the bias can 
fluidize them again, and coarsen their features. This can 
be seen in the transitions from Fig. Illf a) to Fig. Illf b) 
and from Fig.EIc) to Fig.lTTTdL 




FIG. 12: Sedimented interfacial vortex opening out into a 
smoke ring cluster. 3.0 (im spheres in an H = 200 ^im cell, 
(a) A typical vortex ring at 3.4 V. (b) The same cluster at 
2.8 V. 

Increasing the bias still further causes coarsened blobs 
such as those in Fig. ITlT d') to circulate faster and de- 
velop into interfacial toroidal vortex rings, an example of 
which appears in Fig. I12f a~) . These surface- hugging clus- 
ters differ from the free-floating vortex rings of the previ- 
ous section in that they respond irreversibly to changes 
in bias. Reducing the bias on circulating blobs causes 
them to open out into structures resembling smoke rings, 
Fig. H2f bL This morphology can be formed only by first 
increasing and then decreasing the bias, and clearly in- 
herits its topology from the progenitor toroidal cluster. 



8 




FIG. 13: Relaxation behavior of toroidal vortex rings of 3 /im 
silica spheres after a 4 V bias is removed, (a) V: particles 
sedimented onto the lower electrode retain the basic profile of 
their cluster's toroidal geometry. Left in this state, these par- 
ticles would require several minutes to relax to their normal 
diffusive state, (b) Applying a positive bias of 0.8 V levitates 
the particles slightly and allows them to relax to a uniform 
density, (c) Returning to V eliminates any diffusive mo- 
tion and causes the spheres to form sedimented crystalline 
domains, (d) After 3 minutes, the spheres begin to diffuse 
and relax to a homogeneously random state. 



Anomalous dynamics in relaxing colloidal vortices 

Turning off the driving bias after a pattern has formed 
frees the spheres to revert to their equilibrium field-free 
behavior. In most cases, particles sediment into a two- 
dimensional Brownian fluid on the lower electrode. On 
other occasions, and particularly after prolonged runs at 
the upper range of biases, the sedimented particles no 
longer diffuse freely on the lower electrode, but rather 
sit motionlessly where they land. This is the case for 
the particles in Fig. HHf aL which had been organized 
into a colloidal vortex ring before the field was switched 
off. These particles reached the lower electrode rapidly 
enough that their toroidal configuration is still clearly vis- 
ible, and remained unchanged for several minutes. After 
this period, however, the particles began to diffuse nor- 
mally, suggesting that they had not deposited onto the 
lower electrode. This conclusion is supported by the ob- 
servation that a small positive bias suffices to levitate 
the motionless spheres, as shown in Fig. I13f bl . Once the 
positive bias removed, however, the particles sediment 
again to the lower electrode and form small crystals as 
they settle, a typical result appearing in Fig. H3T c). The 
spheres in these self-organized crystals still display no 
Brownian motion for several minutes, after which the re- 
vert to their equilibrium behavior, as shown in Fig. ll3f d). 
This bizarre transient behavior appears not to have been 
described before, and its origin is unknown. 



The islands and labyrinths we have observed somewhat 
resemble the interfacial crystals described in Refs. |l7| 
and 0| , even though the electric field is reversed. 
We conjecture that gravity maintains the particles on 
the lower electrode, while the in-plane electrohydrody- 
namic attraction discovered in Refs. and 0] draws 
them together into islands (Fig. Illf all and labyrinths 
(Fig. UTT c)). Once formed, such close-packed clusters 
would have substantially reduced hydrodynamic drag co- 
efficients j2^, and so would remain sedimented at bi- 
ases that would levitate a single sphere. Applying the 
bias rapidly appears to favor the formation of levitated 
clusters, presumably by disrupting this stabilizing mech- 
anism, while slowly increasing the bias improves the 
chances to form interfacial clusters. 

We have observed free-floating colloidal vortex rings 
coexisting with interfacial labyrinths and islands. This 
demonstrates that the two types of structures can be 
generated under the same conditions, and therefore are 
not distinguished by variations in the experimental pro- 
cedure. 



Dependence on control parameters 

Similar patterns and trends emerge if the bias voltage 
is applied abruptly or increased gradually enough that 
diffuse tumbling clouds form before coalescing into stable 
clusters. The typical number of particles in a steady-state 
cluster tends to increase with the bias. Slightly chang- 
ing the bias once clusters have formed usually does not 
affect this number because merging and splitting clus- 
ters to achieve the new optimal cluster size imposes a 
substantial kinetic barrier. Large bias changes, however, 
can break up all the clusters and reorganize them into 
clusters with the appropriate number of particles. For 
example, groups of three clusters can merge then split 
into cluster pairs with an increase in bias. Reducing the 
applied bias reverses the sequence of transitions, without 
obvious hysteresis. 

The inter-electrode spacing similarly also not to affect 
electrokinetic pattern formation over quite a large range, 
with similar patterns appearing in sample cells at H = 
90, 200, and 390 /im, all with similar crossover voltages. 
Bulk clusters are always observed to float at about h = 
-ff/4, as shown in Fig. El Interfacial patterns form at a 
fraction of this height. Patterns do not form in thinner 
cells, however. For example, we observed that particles 
in a cell H = 40 /im thick were pushed directly to the 
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FIG. 14: Heights of free-floating and sedimented interfacial 
patterns. The error bars represent the variations from run to 
run. In each experiment, all clusters float at the same height, 
with variations well within the error bars. 




FIG. 15: 3.0 /im silica spheres in 0.2 mM NaCl solution. 
H = 200 /im. 4 V bias. Particles repeated clump into dense 
clusters then burst out like fireworks. 

upper electrode and formed no bulk clusters. 

Adding salt slows the particles' motions and favors the 
formation of tumbling clouds instead of quasi-stationary 
clusters. Patterns are blurred in 0.1 mM NaCl solution, 
and disappear altogether in 0.5 mM solution. For ex- 
ample, 3.0 /im silica spheres dispersed in 0.2 mM NaCl 
solution failed to form any quasi-stationary structures, 
but rather entered an irregular cycle of collapsing into 
concentrated domains and then bursting like fireworks 
as shown in Fig. ITSl 

Unlike the previously reported interfacial crystals, 
which form for a wide variety of particle sizes and com- 
positions, the rich panoply of patterns we have de- 
scribed arises appears only for particles meeting specific 



criteria of size and weight. To clarify this, we com- 
pared the electrokinetic pattern formation behavior in 
charge-stabilized dispersions of colloidal silica spheres, 
PS spheres and brominated PS spheres whose densities 
are 1.96, 1.05 and 1.6 g/cm 3 , respectively. The surface 
charge densities of PS and brominated PS are similar and 
do not depend on the pH. Both are more highly charged 
than the silica particles we have discussed so far. Sil- 
ica's charge density, moreover, varies with pH, with an 
isoelectric point at pH 2. 

The behavior of 2.9 /im diameter silica spheres is very 
similar to that of the 3.0 /im diameter spheres we have 
discussed to this point. These also can form various 
bulk and low-floating patterns. Substantially larger silica 
spheres never leave the lower electrode at biases below 
the threshold for bulk electroconvection. Substantially 
smaller, lighter and more diffusive silica spheres create 
less well organized patterns. For example, 1.58 /im sil- 
ica particles usually only form unstable tumbling clouds 
at comparatively low biases before being pushed to the 
upper electrode at higher biases. 

From previous studies jl8| . it is known that electroki- 
netic interactions are weaker for smaller spheres. This 
would tend to reduce the forces binding patterns such 
as colloidal vortex rings together. The smaller spheres' 
more vigorous diffusion also would tend to hinder the 
localized clusters' formation. It also suggests a possible 
role for gravity, which can be highlighted by considering 
spheres comparable in size to the pattern-forming silica 
particles that are composed of different materials. 

Unlike silica spheres, the nearly neutrally-buoyant PS 
spheres fail to form any bulk patterns at any accessi- 
ble bias over the entire diameter range from 0.3 /im to 
7.0 /im. These highly charged spheres should have com- 
parably strong electrokinetic interactions to their silica 
counterparts. Their diffusivity also should be similar. 
Their failure to form bulk patterns demonstrates that 
the gravity plays a central role in the pattern formation 
we have described, a conclusion consistent with the ob- 
servation that no bulk patterns form at negative bias. 

Dense brominated PS spheres 1.6 /im in diameter can 
form tumbling clouds similar to those formed by compa- 
rably sized silica spheres. Whereas clouds of 1.58 /im sil- 
ica spheres tumble vigorously, however, clouds of bromi- 
nated PS particles have no detectable internal convec- 
tion. Such quiescent itinerant clouds also appear for dis- 
persions of 3 /im silica spheres but only at much lower 
concentrations. In both cases, each drifting cloud is 
about 50 /im across and its constituent particles are sep- 
arated by about 10 /im. 

Stable clusters also do not form at very high or very 
low particle concentrations. Dispersions at tj> > 200% 
form chaotic clouds at low bias and jammed clusters 
at high bias. The absence of cyclic motion in low- 
concentration dispersions demonstrates that the flow pat- 
terns we have observed do not arise in the electrolyte 
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alone under the relevant range of conditions, but rather 
require the spheres' presence. The patterns' concentra- 
tion dependence further highlights the role of coopera- 
tive electrohydrodynamic interactions in establishing the 
distinct patterns. Isolated particles' behavior provides 
insights into these forces that help to clarify how such 
patterns form. 



Single-particle behavior 

To best highlight the role of collective behavior in elec- 
trohydrodynamic pattern formation, we studied the sim- 
plest case: a single silica sphere in an electrolysis cell. A 
single charged colloidal sphere of buoyant mass m and 
charge q in a constant electric field not only experiences 
a downward gravitational force mg and an electrostatic 
force qEi but also more complicated electrohydrodynamic 
and osmotic forces due to the electrolyte's response to the 
field. When the lower electrode is negative, the positive 
screening cloud of the sphere flows downward. These 
flowing ions, including OH~ and H30 + created by hy- 
drolysis, entrain a flow of water. The friction between 
this electrohydrodynamic flow and the sphere drags the 
sphere downward. The sphere's cloud of screening ions 
is deformed not only by the electric field, but also by the 
hydronium flux. This downward osmotic force and the 
downward electrohydrodynamic drag are second-order 
responses to the electric field, and therefore should be 
smaller than the electrostatic levitating force qE. 

Straightforward considerations allow us to estimate the 
bias needed to levitate a particle against gravity. A sin- 
gle sphere sediments through water at speed v g — mgb — 
a 2 (p — po)g/(18r]), where po = 1 gm/cm 3 is the den- 
sity of water and rj = 1 cP is its viscosity. For 3 /im 
diameter silica spheres, v g w 5 /im/sec. The same parti- 
cle's electrophoretic velocity saturates at v e « eoS (V ~ 
Vt)/(HKi]), where k _1 is the electrolyte's Debye-Hiickel 
screening length, and eo is the proton charge |l|. If we 
choose the surface charge density T, p h =5 ~ 5000 pm~ 2 
[isl l and kT 1 = 0.1 /im, these two velocity scales would 
be equal at about V — V t = 0.1 V. Thus, considering the 
neglected electrohydrodynamic and osmotic drags, the 
threshold for levitation should be at least 0.1 V higher 
than the water decomposition voltage. This is consistent 
with our experiments. 

In fact, the pH is not uniform across the cell. Fig. 1161 
schematically shows how the profiles of pH, sphere charge 
q, electric field E, and electric force qE vary with height 
in the cell. Even though the electric field is uniform 
in the bulk, the electrokinetic force levitating a charge- 
regulating silica sphere increases with height. When this 
force competes with gravity, therefore, it is possible to 
obtain two equilibrium heights, one near the lower elec- 
trode when gravity dominates, and another near the up- 
per electrode where the electrokinetic force dominates. 



(a) 
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FIG. 16: Gravitational and electrostatic forces on a silica 
sphere, (a) Schematically represents the vertical pH gradi- 
ent established by the dissociation of water during hydroly- 
sis, (b) The charge q on a sphere tracks variations in pH. 
(c) The electric field E is independent of height outside of 
the Debye-Hiickel screening layers, (d) Even so, the charge- 
regulating spheres experience an upward force that increases 
with height in the cell. This force opposes the force due to 
gravity. 



The situation is more complicated for two spheres 
because hydrodynamic forces between the spheres and 
bounding surfaces are long-ranged. Reference |2s| shows 
the numerically calculated streamlines for one or two 
electrophoretically levitated charged spheres. Two par- 
ticles sinking towards the lower electrode create a back 
flow which tends to drive them apart |3f3 |. Conversely, 
two electrically levitated particles will attract each other. 
This electrohydrodynamic interaction at least heuristi- 
cally explains why an initially homogeneous colloidal 
fluid breaks up into discrete clusters as it is driven away 
from the lower electrode. 

We observed that the threshold for levitating a single 
sphere is lower than that for levitating an entire mono- 
layer. This most likely reflects the reduced hydrodynamic 
drag coefficient for many spheres in a monolayer, as com- 
pared to that for a single sphere . The increased hy- 
drodynamic attraction generated by the intensified ionic 
fluxes in the more concentrated dispersions would tend 
to draw these sedimented spheres together into clusters. 
This, however, tends to reduce their overall hydrody- 
namic drag |29j , causing them to rise less than they oth- 
erwise might have. 
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FIG. 17: Height attained by a single 3 /im diameter in a 
H = 200 /im cell as a function of applied bias. The voltage 
is increased and then decreased step by step. Inset: The 
threshold voltage required to levitate spheres. 



most diffuse clouds is about 10 /im and decreases with 
increasing bias. These field-induced attractions are sur- 
prisingly long-ranged, extending to roughly 50 fim. At 
biases above roughly 5 V, even isolated spheres travel in 
circular orbits whose diameters range from 50 to 150 /im 
in the vertical plane. This motion is due to electrocon- 
vection in the underlying electrolyte. It almost certainly 
plays no role in forming the highly organized and sub- 
stantially smaller-scale patterns discussed above. 

Generally speaking, single particles are not stably lev- 
itated into the bulk of the electrolyte, and they certainly 
do not trace out the complex trajectories characterizing 
microscopic patterns such as those in Figs. Qj]- El These 
observations help to confirm that such patterns do not 
arise from electroconvection in the electrolyte alone. We 
are left to conclude that the spheres play an active role 
in dynamical pattern selection. 



DISCUSSION 



To compare our system's behavior with these predic- 
tions, we performed experiments on 3 /im diameter sil- 
ica spheres from a monolayer with areal density below 
10~ 4 /im -2 , i.e. fewer than three particles per frame. 
Abruptly applying a bias causes these isolated spheres to 
jump off the lower electrode during the transient in which 
vertical ionic concentration gradients are established. Be- 
low a threshold voltage, the particles eventually return 
to the lower wall. Above this threshold, they rise to the 
top wall and remain suspended. Fig. 1171 shows how the 
single sphere's height changes with applied voltage as the 
voltage is first increased and then decreased. The point 
at which a single sphere is electrolevitated to the up- 
per wall is hysteretic: levitated particles only return to 
the lower wall when the voltage is lowered substantially. 
The explanation for this hysteresis also can be found in 
Fig- El At the critical bias, |gE(z)| > mg for all z in 
the bulk. By contrast, a particle can only fall from the 
upper electrode, when the field falls below threshold at 
the upper electrode \q~E(H)\ < mg. Therefore, the two 
thresholds are different. The inset to Fig. 1171 shows this 
hysteresis more clearly at three different salt concentra- 
tions. Such hysteresis supports the contention that the 
levitating force acting on silica spheres increases with the 
height h. 

For salt concentration above 1 mM, a single particle 
can never be forced to the upper electrode. This is con- 
sistent with the observation that particles move less vig- 
orously and cooperative structures are harder to form at 
high ionic concentrations. 

Attractions between particles are detectable only at 
relatively low particle concentrations (e.g. <fi ~ 10~ 4 ) 
and high applied voltages (e.g. 2 V). In such cases, par- 
ticles attract each other as they levitate and so form tum- 
bling clouds. The typical inter-particle separation in the 



We have demonstrated that many distinctive patterns 
reproducibly emerge from the interplay of gravity and 
electrohydrodynamic forces in charge-stabilized colloid. 
At biases above the threshold for hydrolysis but below 
the onset of electroconvection, only fluxes of HaO + and 
OH~ ions flow through the system. Electroneutrality 
favors overlap of these fluxes; however, "friction" be- 
tween opposite flows is minimized if these flows are sep- 
arated. The optimal balance between these effects sep- 
arates the fluxes by a characteristic lengthscale compa- 
rable to the Debye-Hiickel screening length, or roughly 
0.1 /im. This is much smaller than the extent of the ex- 
perimentally observed clusters whose diameters extend 
to roughly 100 /im. Electroneutrality suppresses the for- 
mation of such large-scale convective structures in the 
absence of the charged spheres. In other words, are not 
simply passive tracers for a pre-existing flow in the elec- 
trolyte, but rather they and their charges are essential 
for establishing the convective flows at the heart of the 
observed patterns. 

There are three general features about the dynamics of 
all clusters. (1) They have the same circulating direction 
as shown in Fig. (2) The circulation speed increases 
with the bias as clusters tend to become more compact. 
(3) The number of particles in each cluster increases with 
the bias. The bias-induced attraction of the second fea- 
ture is crucial to understand the clusters' formation. The 
breathing mode in Figs. 151 and !15l mav result from compe- 
tition between this as-yet unexplained attraction and a 
combination of screened-Coulomb repulsions and induced 
dipolc repulsions among spheres. 

While the sequence of dynamical transitions from sedi- 
mented layer to tumbling clouds to highly organized clus- 
ters is reproducible, the crossover potentials depend sen- 
sitively on the changing electrochemical properties of the 
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electrodes, and therefore can vary by as much as ±0.4 V 
from run to run. 

SUMMARY 

We have reported a large family of previously unrecog- 
nized self-organized colloidal patterns that form in con- 
stant vertical electric fields at biases just above the de- 
composition voltage of water. Our simple electrolysis cell 
gives rise to a complicated coupled system of hydrody- 
namic flows and ionic fluxes that are shepherded by the 
very colloidal particles they transport. We have discussed 
the electrohydrodynamic forces on a single particle in this 
system, and suggested how many-body coupling might 
give rise to the various dynamical structures we have 
observed. At low voltages, colloidal spheres in a spe- 
cific range of charges, densities and mobilities are found 
to form various microscopic quasi-steady-state clusters 
through the competition of gravity and electrohydrody- 
namic levitating forces. These patterns falls into two 
categories: clusters levitated into the bulk and clusters 
localized near the lower electrode. Their structures and 
evolutions are different, but all share some general fea- 
tures, such as persistent circulation and coarsening with 
increasing voltage. The simple behavior of a single parti- 
cle in the field demonstrates that these microscopic clus- 
ters are formed cooperatively. Patterns' dependence on 
such control parameters as the composition and number 
density of particles, cell thickness H , salt concentration 
and applied voltage range have been discussed. 
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Petroleum Research Fund of the American Chemical So- 
ciety. 
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